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Abstract 

Background: Microsatellites or simple sequence repeats (SSRs) are highly informative molecular markers for various 
biological studies in plants. In spruce (Picea) and other conifers, the development of single-copy polymorphic genomic 
microsatellite markers is quite difficult, owing primarily to the large genome size and predominance of repetitive DNA 
sequences throughout the genome. We have developed highly informative single-locus genomic microsatellite markers in 
black spruce (Picea mariana) and red spruce (Picea rubens) using a simple but efficient method based on a combination of 
AFLP and microsatellite technologies. 

Principal Findings: A microsatellite-enriched library was constructed from genomic AFLP DNA fragments of black spruce. 
Sequencing of the 108 putative SSR-containing clones provided 94 unique sequences with microsatellites. Twenty-two of 
the designed 34 primer pairs yielded scorable amplicons, with single-locus patterns. Fourteen of these microsatellite 
markers were characterized in 30 black spruce and 30 red spruce individuals drawn from many populations. The number of 
alleles at a polymorphic locus ranged from 2 to 18, with a mean of 9.3 in black spruce, and from 3 to 1 5, with a mean of 6.2 
alleles in red spruce. The polymorphic information content or expected heterozygosity ranged from 0.340 to 0.909 
(mean = 0.67) in black spruce and from 0.161 to 0.851 {mean = 0.62) in red spruce. Ten SSR markers showing inter-parental 
polymorphism inherited in a single-locus IVlendelian mode, with two cases of distorted segregation. Primer pairs for almost 
all polymorphic SSR loci resolved microsatellites of comparable size in Picea glauca, P. engeimannii, P. sitclnensis, and P. abies. 

Significance: J\r\e AFLP-based microsatellite-enriched library appears to be a rapid, cost-effective approach for isolating and 
developing single-locus informative genomic microsatellite markers in black spruce. The markers developed should be 
useful in black spruce, red spruce and other Picea species for various genetics, genomics, breeding, forensics, conservation 
studies and applications. 
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Introduction 

Microsatellites or simple sequence repeats (SSR) are highly 
informative genetic markers. Due to their co-dominant inheri- 
tance, high polymorphism, reproducibility and transferability 
across related species, microsatellites have been widely used for 
various genetics, genomics, breeding, genetic resource conserva- 
tion, and forensics studies and applications, e.g., [1-8]. Although, 
recendy with the use of next-generation sequencing (NGS) 
technologies, genome-wide single nucleotide polymorphism 



(SNP) and other markers have become available, microsatellite 
markers are stiU highly suitable for population and conservation 
genetics studies and forensics applications, especially for deter- 
mining neutral population genetic processes and dynamics. In fact, 
microsatellites were demonstrated to provide more precise 
information on population genetic structure than SNPs in four 
diflFerent taxa [9] . Also, highly polymorphic co-dominant markers, 
such as microsatellites, are needed to align genetic maps. 
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Various methods and strategies have been used to isolate 
microsatellite-containing sequences in plants and animals, which 
are primarily based on the development of microsateUite-enriched 

and non-enriched genomic libraries [10]. Also, availability of a 
large number of expressed sequence tag (EST) sequences in the 
database made it possible to identify microsatellite sequences and 
develop cDNA-based microsatellite markers in agricultural and 
horticultural plants, e.g., [11-13] and conifer trees [14,15], The 
advantages associated with EST-derived microsatellite markers 
include high probability of such markers being single-locus types 
and their ability to detect polymorphism directiy in expressed 
genes. However, EST-derived microsatellite markers show lower 
polymorphism than genomic microsateUites [7]. This is possibly 
due to lower mutation rates, conservation and slow evolution of 
gene sequences and selection constraints on them. Therefore, 
information obtained from EST-derived microsateUites is likely to 
be different from that obtained from faster-evolving genomic 
microsateUites. Also, the expressed portion of the genome is quite 
smaU, especially in plants with large genome size; thus, cDNA- 
based microsatellite markers may have limited and non-random 
representation of the genome. Hence, it is important to develop 
genomic microsateUites. More recentiy, massive parallel sequenc- 
ing of the whole transcriptomes and genomes can reveal a large 
number of sequences with microsatellite repeats. 

Due primarily to huge genome size (2-4x10'" base pairs) and 
highly repetitive genomes of conifers [16,17], development of 
highly informative single-locus genomic microsatellite markers 
remains a chaUenge in these long-lived plants. Although a large 
number of microsatellite-containing sequences can now be 
identified by whok; genome- sequc-ncing in conifers [18,19] using 
NGS platforms, the challenge still remains to isolate and develop 
genomic microsatellite markers that show simple single-loi:us DNA 
variant patterns. Genomic microsateUite markers have been 
developed for a number of spruce (Picea) and other conifer 
species using non-enriched and microsateUite-enriched [2,20-23] 
and BAC/YAC genomic libraries [24,25]. However, the number 
of genomic microsateUites showing single-copy single-locus 
patterns remains small because a large proportion of genomic 
microsateUites have shown complex multilocus patterns with a 
frequent presence of nuU alleles. The difficulties and chaUenges 
have been mainly attributed to the large genome size and the 
presence of high proportion of redundant repetitive DNA 
sequences throughout the genomes in conifers [16,17]. In order 
to minimize the problems, attempts were made to isolate 
microsateUite markers from unique DNA sequences, using two 
approaches: low copy enrichment method [26] and under- 
methylated (UM) DNA method following methylation restriction 
enzyme (McrBC) [27]. Of the above two approaches, the l<n\--c()py 
enrichment method was found to be more successful. However, 
both of these methods are stiU time consuming and expensive and 
yielded only a smaU number of highly informative single-locus 
microsatellite markers. Development of informative single-locus 
microsatc-Uitc markers \\ as h()\vc\'("r more successful from cDNA 
or EST sequences in conifers [13-15]. 

Here we report a simple and efficient approach for the 
development of highly informative single-locus genomic microsat- 
eUite markers in two widely distributed and economically and 
ecologically important conifers, black spruce {Picea mariana) and 
red spruce [Picea ruhens). Microsatellite sequences were identified 
from microsateUite-enriched AFLP DNA fragments in black 
spruce. We have determined informativeness of these markers in 
black spruce and red spruce, their Mendehan single-locus 
inheritance and linkage in progeny of two controUed crosses of 



black spruce, and their cross-species transferabUity to four other 
spruce species. 

Materials and Methods 

Plant Material and DNA Extraction 

One of the parents (accession # 59) of a three-generation 
outbred mapping pedigree of black spruce was used for the 
construction of microsatellite-enriched library. Thirty individuals 
each of black spruce and red spruce drawn from 1 2 populations 
were used to determine the informativeness of the microsatellite 
markers. Eight individuals each of white spruce (P. glauca), Sitka 
spruce [P. sitchensis), Engelmann spruce [P. engelmannii), and 
Norway spruce (P. ahies) were used for determining the cross- 
species transferabUity of the microsatellite markers, with the 
exception of one marker RPMSA07 for which five individuals of 
each species were used. The black spruce samples were previously 
obtained from nine populations located at three sites in Manitoba, 
Canada [28] as follows: Pine Falls - natural mature (1), post- 
harvest natural young regeneration (1), and plantation (4); Bissett - 
natural mature (3), post-harvest natural young regeneration (2), 
and plantation (3); and Snow Lake - natural mature (3), post- 
harvest natural young regeneration (6), and plantation (7). The red 
spruce samples were from three provenances included in a range- 
wide provenance test established at the Acadia Forest Experiment 
Station, Fredericton, New Brunswick: Great Smoky Mountain 
National Park, Tennessee-North Carolina boundary (6), Mono- 
ngahela National Forest, Pocahontas County, West Virginia (12), 
and Logan State Forc'st, Pennsylvania (12). The Norway spruce 
samples were ohtainc'd from the Nova Scotia Department of 
Natural Resources, Sitka spruce samples from the British 
Columbia Ministry of Natural Resources, and white spruce and 
Engelmann spruce samples were from Alberta, Canada. The 
parents and 1 10, and 30 progeny of two controUed crosses of 
black .spruce 32x40, and 46x14, respectively, were used to 
examine the inheritance and linkage of microsatellite markers. 
The detaUs on controUed crossing experiments and pedigrees are 
provided in Kang et al. [6]. 

Total genomic DNA was extracted from the needle tissues of 
each individual using Qiagen DNeasy plant mini kit according to 
the manufacturer's instructions (Qiagen, Toronto, ON, Canada). 
The quality and quantity of DNA was evaluated using a 
spectrophotometer and agarose gel electrophoresis. 

Isolation of Microsatellite DNA Sequences and 
Development of SSR Markers 

An AFLP-based microsatellite-enriched Hbrary was constructed 
using total genomic DNA of black spruce accession 59 and 
applying the magnet beads method [29,30], with minor modifi- 
cations. Genomic DNA (~600 ng) was digested with two different 
restriction endonucleases iicoRI and Msel at 37°C for 2 h and 
ligated with £coRI and Msel adapters using the AFLP Analysis 
System I (Invitrogen Life Technologies, Burlington, ON, Canada). 
An initial PCR reaction was conducted in a total volume of 200 jil, 
containing 2.5 mM MgCl,, 0.16 mM each of dNTPs, 10 pmoles 
of £coRI adaptor-specific primer (5'-GACTGCGTACCAATTC- 
3') and 60 pmoles of Msel adaptor-specific primer (5'-GAT- 
GAGTCCTGAGTAA-3'), 10 unit of Taq polymerase (Fermen- 
tas. Fisher Scientific, Ottawa, ON, Canada) and 10 |.ll of purified 
adaptor-ligated DNA as templates. Amplification was performed 
with 10 cycles each of denaturation at 90°C for 30 sec, annealing 
at 56°C for 60 sec and extension at 72°C for 60 sec. The 
amplification products were hybridized with 100 pmoles of 3'- 
biotin-labeUed (AG) 15 oligonucleotides, and the target DNA 
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fragments were recovered using Dynabeads M-280 Streptavidin 
(Dynal Biotech, Life technologies, Burlington, ON, Canada) in a 
STEX buffer (100 mM NaCl, 10 mM Tris-HCl, 1 mM EDTA 
and 0.1% Triton- 100, pH 8.0) for 30 min [29,30]. We used 
(AG) 15 probe for microsatellite enrichment because AG repeats 
were found to be most abundant in white spruce [2] . The enriched 
DNA fragments were amplified again with 25 cycles using the 
same cycling conditions. The purified PGR products were inserted 
into the pCR2.1 TOPO vectors and then transformed into One 
Shot Competent Cells provided in the TA cloning kit (Invitrogen 
Life Technologies, Burlington, ON, Canada), following the 
manufacturer's instructions. The bacteria were plated on LB agar 
plates containing X-gal and 100 |Xg/ml ampiciUin. The positive 
clones (white and some light blue colonies) were screened by PGR 
using the Universal Ml 3 Reverse Primer and the T7 Promoter 
Primers (Invitrogen Life Technologies, Burlington, ON, Canada). 
Colony blot hybridization was carried out using biotenylated 
(AG) 15 probe and BrightStar BioDetect non-isotopic detection kit 
according to the user's manual (Ambion, Life Technologies, 
Burlington, ON, Canada) to further screen and identify the 
plasmids containing microsatellite sequences. 

The putative microsatellite-containing clones were grown 
overnight in 3 ml of LB broth with 50 |J,g/ml ampiciUin. 
Subsequently, plasmids were extracted using a QIAprep Spin 
Miniprep Kit (Qiagen, Toronto, ON, Canada). The purified DNA 
fragments were sequenced using the Universal M13 Reverse 
Primer labelled with IRdye 800 and the T7 Promoter Primer 
labeled with IRDye 700 with Thermo Sequenase fluorescent- 
labeled primer cycle sequencing kit (Amersham Pharmacia 
Biotech, Quebec, Canada). 

The sequence data were ansdyzed using the E-SEQ_v2 program 
and the Align IR program (LICOR, Lincoln, Nebraska, USA). 
Primers were designed for the non-redundant SSR-containing 
sequences using Primer 3 program [31]. The oligonucleotide 
primers complementary to the regions flanking the identified SSR 
repeat motifs were synthesized incorporating IRDye 700 or 800 
labels by the Integrated DNA Technologies (CoralviUe, lA, USA), 
and used for amplification by PGR. 

PCR Amplification, Optimization and Visualization of 
Microsatellites 

PGR amplification for all microsatellite markers was performed 
in a 20 nl reaction volume, with 100 ng of purified spruce 
genomic DNA, 1.5 mM of MgClj, 200 |J,M of each of dNTPs, 
200 ng of BSA, 1 U of Taq DNA polymerase (Fermentas, Fisher 
Scientific, ON, Canada) and 7.5 pmole of each forward and 
reverse primers. To optimize the amplification conditions for each 
microsatellite locus, the following temperature cycling parameters 
were tested: denaturation for 3 min at 94°C, followed by two 
cycles of 30 sec each at 94°G, 55-60°C and 72°C; and 38 cycles 
each of 15 sec at 94°C, 55-60°G and 72°C, with a final extension 
step of 10 min at 72°C. For the primers yielding no, weak or 
multi-locus amplification products using the above PCR profiles, a 
touchdown protocol was used: initial denaturation at 94°G for 
3 min followed by two cycles each of 30 sec at 94°C, 60°G and 
72°C; followed by 13 cycles each of 94°C for 15 sec, 60 i 54°C for 
15 sec (0.5°C per cycle) and 72°C for 15 sec; and 30 cycles each of 
15 sec at 94°G, 54°C and 72°C, with a final cycle at 72°C for 
10 min. 

After amplification, the PGR products were separated on a 
6.5% polyacrylamide gel at 1500 volts for 2.5 h on the LICOR 
sequencing system IR 4200 (LiCor, Lincoln, NE, USA) with 
IRDye 700 and 800 labelled size standards (IR 50-350 bp). The 
genotypes of spruce individuals were determined by their aUeUc 



constitution manually. The DNA fragment size of individual 
alleles was determined based on the size standard used. 

Informativeness, Inheritance and Linkage of 
Microsatellite Markers 

The informativeness of the microsatellite loci was evaluated in 
black spruce and red spruce by determining the number of alleles 
and their frequency, observed heterozygosity(H„), and expected 
heterozygosity (HJ. Also polymorphism information content (PIG) 
for each locus was calculated as follows [32]: 1 — Yl^^ where P; 
is the frequency of the ith allele in 30 individuals of a species 
examined. The PIG is the same as expected heterozygosity (H,.). 

The inheritance, segregation patterns and linkage of the 
polymorphic microsatellite markers were examined in two three- 
generation outbred pedigrees, which consisted of the parents and 
F2 progenies of two controlled crosses (32 x40 (1 10 progeny), and 
46x14 (30 progeny) of black spruce. The linkage relationships 
among the seven microsatellite loci (RPMSA04, RPMSA07, 
RPMSA09, RPMSA12, RPMSA13, RPMSA17, and RPMSA33), 
showing interparental polymorphisms in the 32x40 cross have 
been previously determined [6]. In the 46x14 cross, the two-point 
Hnkage was determined among five informative loci (RPMSA04, 
RPMSAll, RPMSA12, RPMSA26, and RPMS33), and between 
RPMSA13 and RPMSA22 by investigating joint segregation and 
independent assortment of their microsatellite DNA variants. Due 
to certain logistics reasons, three different sets of 30 F-i progeny of 
the 46x14 cross were used for determining the inheritance and 
linkage of the microsatellite markers: set 1 (RPMSA04, 
RPMSAll, RPMSA12, RPMSA26, and RPMS33), set 2 
(RPMSA07), and set 3 (RPMSA13, and RPMSA22). Therefore, 
the pair-wise finkage relationships among all eight markers 
showing interparental polymorphism in this cross could not be 
tested. The null hypothesis that the two loci in a pair are 
independentiy assorted, was tested by the goodness-of-fit test. 
Since the expected number of progeny in genotypic classes for five 
microsateUite pairs (RPMSA04-RPMSA 1 1 , RPMSAll- 
RPMSA12, RPMSAll -RPMSA26, RPMSAll- RPMSA33, and 
RPMSA13-RPMSA22) was less than the minimum of 5 required 
for a vaUd test [33], we pooled adjacent two or four progeny 
genotypic classes to obtain a joint class with an expected progeny 
number >5 to perform the test [33]. This resulted in a total of 
two or four joint progeny classes. 

Cross-species Transferability of Microsatellite Markers 

The number of alleles observed at each locus in eight individuals 
(five individuals for RPMSA07) of Norway spruce, Sitka spruce, 
white spruce, and Engeknann spruce were determined to test the 
cross-species transferability of the microsatellite markers developed 
from black spruce. 

Results 

Construction of Microsatellite-enriched Library and 
Identification of SSR-containing Sequences 

The microsatellite-enriched library was constructed from 
genomic DNA of black spruce after restriction with a rare and a 
frequent endonucleases and enrichment with biotinylated (AG) 15 
oligonucleotides. One-thousand twenty-nine white and light blue 
colonies were isolated. Of these, 311 were identified as positive 
potentially containing DNA inserts after the first colony PGR 
screening with the T7 Promoter and the Universal Ml 3 reverse 
primers. The second screening with the (AG)i5 oligonucleotides as 
probes using non-isotopic detection method resulted in excluding 
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Table 1. Polymorphic microsatellite markers/loci showing simple single-locus patters; and their repeat motifs, primer sequences, 
annealing temperatures, and GenBank accession numbers. 





Microsatellite locus 


Repeat type 


Primers (5 -3 ) 


Annealing temperature (°C) 


GenBank Accession No. 


RPMSAOl 


(TA)9(GA)7(GATA),3 


GGAGCTAAACACATTTGGTACAGG 


60 i 54 


KJ847201 


GAAACCATTGATGTGGGTTG 


RPMSA04 


(CT)23 


GTGAGATTTTGGCAGCAACA 


60 


KJ847204 


TGATCACCCTTGCTCAAAGA 


RPMSA05 


(GA)2o 


CCCTATTCCCACTTGAAATCC 


60 


KJ847205 


CTTATGGGCTCCACCACACT 


RPMSA06 


(CT)2(GA)6 


CTCGACAGACCCCTCTTTTG 


60 J, 54 


KJ847206 


CAAGTCTCGGTTCTCCTCCA 


RPMSA07 


(CT),6 


TGAAGAAGCAAGTGGGCTCT 


60 i 54 


KJ847207 


TGCATTGATCTCTCCCCTTT 


RPMSA09 


(AG)a 


CACCTCAGTTCACACCTGCT 


60 


KJ847209 


TTCCTCTCCCAAGAATGTGC 


RPMSAU 


(GA)2, 


GACCCTAGATTTTGGGGTAT 


60 


KJ847211 


CCCCCTCTCAGTAATCCAAC 


RPMSA12 


(TC)5CC(TC)4 


AACGAGGTTCATCCCATCTG 


60 i 54 


KJ847212 


TACGCTCAATGTCGATGAGG 


RPMSA13 


(GA)9 


AACCATGAAACCCTAGCGACT 


60 


KJ847213 


TGAGGACTTAGGCCCACATT 


RPMSA15 


(GA)7(AGAAT)(GA)4 


ATCGATAGGCTTGCAAGAGG 


60 i 54 


KJ847215 


TGTGCCCTCATGTGCTATGT 


RPMSA17 


(CT),2 


CAACGACTGCAACTGGGTACT 


60 


KJ847217 


CAACCATAGACACGCAACCA 


RPMSA19 


(CT),2 


TAGCCAATACAATGCCAAGG 


60 


KJ847219 


ATCAGAGCGAAGTTTGGAG 


RPMSA22 


(TC)6/(TC),4/(GA)4 


GCACGTGCATGTTCTCTGTC 


58 


KJ847222 


TGCATGCAGATGAATGAGAG 


RPMSA26 


(CT)2C(CT)9 


GCTGTAGGGTTGATATTTGC 


65 I 60 


KJ847226 


TGTGTGAGAGATAAGTGTTGAG 


RPMSA27 


(TC)22(TA),5 


ATATTCGAATGAGAGAGCAATC 


55 


KJ847227 


TGTAGGCCCATGATAATGTA 


RPMS/lii 


(GA), 


ACACACATGAACACATGAGC 


60 i 55 


KJ847233 


GCTGTATGGATTCCGTATGA 



doi:l 0.1 371 /journal.pone.Ol 03789.t001 



203 clones from the putative SSR-containing clones. Sequencing 
of the remaining 108 putative positive clones revealed 94 unique 
sequences containing microsatellite repeats with mono- (8.5%), di- 
(87.3%), tri- (3.2%), and tetranucleotide (1%) motifs (Table SI). 

The DNA insert sizes of the 94 microsatellite-containing clones 
ranged from 148 bp to 1,273 bp, with an average of 462 bp. The 
majority of clones (77.4%) contained a moderate insert size of 
200-600 bp, while five clones (BS-013, 019, 026, 236 and 245) 
contained larger DNA inserts (>800 bp) and three clones (BS-039, 
1 19 and 241) smaller inserts (<200 bp). Thirty-nine (41.5%) of the 
identified microsateUites were perfect, 27 (28.7%) imperfect and 28 
(29.8%) compound (Table SI). In terms of repeats, GA/CT motif 
was the most common, accounting for 84.6% of the perfect 
repeats. One clone (BS-236) contained a simple AT motif and the 
other two TG/CA motifs. A perfect trinucleotide (CTT) repeat 
motif was also found in one sequence (BS-120). Among the 
compound microsateUites, 25 were composed of two repeat types 



and three of three repeat types. Twenty-one out of 28 compound 
microsateUites were interrupted with one or more nucleotides 
(Table SI). The overall repeat motif number ranged from 4 to 54 
units (Table SI). 

Microsatellite Markers and Their Informativeness in Black 
Spruce and Red Spruce 

Primer pairs were designed from the flanking regions of 
microsateUites for selected 34 of the 94 weU-characterized 
microsatellite-containing sequences (Table S2). Each of these 
primer pairs were initially tested for amplification of microsateUite 
DNA variants in three black spruce genotypes (59, 63 and N2). A 
variety of PGR conditions were tested. Of this set, 22 (64.70%) 
primer pairs yielded single-locus microsatellite variant patterns. 
The remaining 1 3 SSR primers either did not yield any amplicon 
or showed multiple DNA fragment patterns. Among the 22 
microsatellite primer pairs resolving microsateUites with simple 
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single-locus patterns, 16 primer pairs produced consistent and 
clear scorable amplicons in 30 black spruce individuals and 15 
primer pairs in 30 red spruce individuals (Table 1; Table 2; 
Fig. 1). The optimal annealing temperatures for each of the 16 
primer pairs are listed in Table 1 . 

In black spruce individuals, two loci, RPMSAOl and 
RPMSA06, both compound microsatellites (Table 1) were found 
to be monomorphic (Table 2). The numbers of alleles detected at 
the polymorphic markers were quite variable ranging from 2 
(RPMSA26) to 18 (RPMSA27), widi an average of 9.3 aUeles per 
locus (Table 2). The DNA fragment sizes of alleles ranged from 
105 bp to 267 bp (Table 2). Among the red spruce individuals, the 
RPMSAO 1 locus was monomorphic, while the RPMS A05 primers 
produced multilocus patterns. The number of aUeles ranged from 
3 to 15 at a locus in red spruce (Table 2), with a mean of 6.2 aUeles 
per locus. The most informative locus in red spruce in terms of the 
number of alleles was RPMSA27, with 15 alleles detected in the 30 
individuals. The frequencies of individual alleles at the polymor- 
phic loci in 30 individuals of black spruce are provided in Table 
S3, whereas that for red spruce in Table S4. Most of the alleles had 
low frequencies as expected for highly polymorphic loci. In black 
spruce, the allele frequencies at a polymorphic locus ranged from 
0.017 at 13 loci to 0.80 at RPMSA26 (Table S3). In red spruce, 
frequencies of individual alleles ranged from 0.017 at seven 
polymorphic loci to 0.588 at RPMSA09 (Table S4). Red spruce 
individuals had 30 alleles not detected in black spruce samples, 
whereas black spruce samples had 70 alleles not detected in red 
spruce samples at 15 loci (RPMSA05 was excluded because its 
primers produced multilocus patterns in red spruce and genotypes 
of individual samples could not be scored) (Table S3; Table S4). 
The expected heterozygosity or PIC values for the polymorphic 
microsatellite loci varied from 0.34 (RPMSA26) to 0.909 
(RPMSA05), with an average of 0.67 in black spruce, and from 
0.161 (RPMSA04) to 0.851 (RPMSA13), with an average of 0.62 



in red spruce (Table 2). The observed heterozygosity values for the 
polymorphic microsatellite loci varied from 0.100 (RPMSA12) to 
0.900 (RPMSAll and RPMSA15), with an average of 0.401 in 
black spruce, and from 0.133 (RPMSA22) to 0.900 (RPMSA15), 
with an average of 0.381 in red spruce (Table 2). 

Inheritance and Linkage of SSRs in Black Spruce 

Ten microsatellite loci were found to be polymorphic between 
the parents of one or both controlled crosses (Table 3; Fig. 2). At 
each of the 1 0 loci, the progenies of each cross segregated into two 
to four genotypic classes as expected for single-locus Mendelian 
inheritance patterns (Table 3). However, significant deviations 
from the expected MendeKan progeny ratios were observed in two 
cases. Segregation of alleles at RPMSA04 in one of the two 
controlled crosses (46 x 1 4) conformed to the expected Mendelian 
ratios. However, the observed ratio in the progenies of the other 
cross (32x40) departed from the Mendelian expectations as 
evaluated by x^-tests, PsO.05 (Table 3). The second case was with 
RPMSA17 with four segregating alleles (Table 3). 

In the 32x40 pedigree, seven microsatellite loci showed 
independent assortment, and six of them were previously mapped 
on to six linkage groups (LG) [6]: LGOl: RPMSA13; LG02: 
RPMSA09; LG04: RPMSA12; LG05: RPMSA33; LG06: 
RPMSA17; LG 07:RPMSA4. In the 46x14 cross, die Chi-square 
goodness-of-fit tests indicated that the observed numbers of 
genotypes were in agreement with those expected for independent 
assortment of two loci in each of the 1 1 two-locus combinations 
tested (Table S5). 

Cross-species Transferability of Black Spruce 
iVlicrosatellite Markers in other Picea Species 

Reproducible amplifications of the microsateUites of most of the 
15 loci were obtained in Norway spruce, Sitka spruce, Engeknann 
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E2 PI- 



PROGENY 



RPMSA04 




Figure 2. Inheritance of microsatellite DNA variants in black spruce. Inheritance and segregation pattern of microsatellite markers at (a) 
RPI\/ISA04, and (b) RPMSA11 in F2 progeny of 46x14 controlled cross. The genotypes of the parents and progeny are provided in Table 3. 
doi:1 0.1 371/journal.pone.01 03789.g002 



spruce, and white spruce (Table 4). The fragments amplified from 
the four species were in the size range as obser\'ed in black spruce. 

Among the 15 primers tested, nine yielded polymorphic 
amplicons in at least two species (Table 4). The RPMSA33 
primers yielded polymorphic PGR products in Engelmann spruce 
and white spruce, a monomorphic product in Sitka spruce and no 
product in Norway spruce. The RPMSA0,5 primers detected 
several loci in Norway spruce, Sitka spruce, Engelmann spruce, 
white spruce and red spruce. The RPMSAOl SSR locus was 
monomorphic in all spruce species tested, including black spruce 
and red spruce. Interestingly, RPMSA06, RPMSA9, and 
RPMSA19 resolved polymorphic microsatellite DNA variants in 
black spruce and/ or red spruce, but these loci were monomorphic 
in the individuals of other four species. While two or three alleles 
were detected at RPMSA26 in black spruce and red spruce 
samples (Table 2), the primers for this locus did not show any 
amplification in the four species. 



Discussion 

Isolation and Development of Microsatellite DNA IVlarkers 

We have developed highly informative single-locus genomic 
microsatellite DNA markers in black spruce and red spruce using 
an AFLP-based microsatellite enrichment and isolation method. 
This is the first report for microsatellite marker development in red 
spruce. A few genomic microsatellite markers have recendy been 
reported for black spruce [34]. We demonstrate that 14 of the 16 
microsatellite markers characterized are informative in black 
spruce and red spruce. Of these, seven (RPMSAO.^i, RPMA07, 
RPMSAll, RPMSA13, RPMA17, RPMSA27, and RPMSA33) 
were the most informative loci in black spruce, each with over 1 0 
alleles and expected heterozygosity/PIC close to 0.90. In red 
spruce, four loci (RPMAll, RPMSA13, RPMA17, and 
RPMSA27), with 10 or more alleles at a locus were the most 
informative. The levels of allelic polymorphisms of the microsat- 
ellite loci were positively correlated with the repeat numbers 
present in the original clones isolated from black spruce (r = 0.484, 
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Table 4. Transferability of microsatellite markers to other four spruce species: number of alleles detected in selected individuals of 
each species. 



Locus 


Number of alleles 








Norway spruce 


Sitka spruce 


Engelmann spruce 


White spruce 


RPMSAOl 111 1 


RPMSA04 


2 


8 


3 


1 


RPMSA05 


Multilocus 


Multilocus 


Multilocus 


Multilocus 


RPMSA06 


1 


1 


1 


1 


RPMSA07* 


6 


5 


4 


5 


RPMSA09 


1 


1 


1 


1 


RPMSAll 


3 


10 


10 


9 


RPMSA12 


5 


2 


2 


3 


RPMSA13 


6 


9 


9 


7 


RPMSA17 


10 


7 


8 


6 


RPMSA19 111 1 


RPMSA22 


7 


5 


7 


4 


RPMSA26 . _ . 


RPMSA27 


6 


6 


6 


8 


RPMSA33 




1 


5 


4 


Note: 



*Five individuals each of P. abies, P. sitchensis, P. glauca and P. engelmannii were used to amplify the SSR locus RPMSA07, whereas eight individuals per species were 

used to amplify other 14 SSR loci. 

- SSR locus not detected. 

doi:1 0.1 371 /journal.pone.Ol 03789.t004 



^ = 0.04) and red spruce (r = 0.480, ^ = 0.04). These re.sults are 
consistent with the positive correlation observed between the 
repeat length and the polymorphism in white spruce [2] , Norway 
spruce [22], and other plants [35,36]. In contrast, no such positive 
correlation between the microsatellite repeat length and resultant 
polymorphism was observed in some plant species [37,38]. 
Microsatellite loci with compound repeats are generally believed 
to show lower polymorphism than loci with simple repeats. While 
this was true for two (RPMSAOl, and RPMSA06) of the six 
microsatellite loci with compound repeats in our study, the locus 
RPMSA27 with compound CT and AT repeats showed the 
highest polymorphism in black spruce and red spruce (Table 2). 

Several approaches have been used to isolate and develop 
microsatellite DNA markers from spruce species (Table 5). The 
comparison with other approaches (Table 5) suggests that the 
AFLP-based approach that we have reported here, provided the 
most efficient method for isolating and developing highly 
informative single-locus microsatellites with simple patterns in 
spruce. The success rate was over 250% higher than that was 
observed for microsatellite marker development in black spruce 
from standard microsatellite-enriched genomic libraries in our 
own lab (Table 5; Joy and Rajora, unpubhshed). It is worth noting 
that the success rate was even higher than that observed for 
isolating and developing microsatellites from cDNA and EST 
sequences [15] or from low-copy sequences selected from 
microsatellite-enriched genomic libraries [39] in spruce. We used 
£coRI as one of the two restriction enzymes, which is methylation 
sensitive. The repetitive elements in the genome are highly 
methylated [40]. Therefore, methylated sites in the genome may 
have been preferentially excluded by using £coRI and this may 
have enriched single- or low copy sequences in the microsatellite- 
enriched AFLP fragments library. Similar AFLP-based approach 



has also been successfully used to develop microsatellite markers in 
pear [Pyrus pyrifolia) [29] and peach [Prunus persica) [30]. 

An AFLP-based method named Fast Isolation by AFLP of 
Sequences Containing Repeats (FIASCO) was recently introduced 
[10], which was believed to be a simple, rapid and effective 
method for microsatellite isolation from genomic DNA de-novo. 
However, this method did not prove effective in conifers as the 
percentage of microsatellite-containing clones was relatively low 
(between 29.5 to 43.1%) [23,41,42]. The FIASCO method does 
not apply a microsatellite enrichment step, whereas our method 
employs microsatellite enrichment followed by double screening 
for microsatellite-enriched clones. In our study, the enriched 
genomic library was screened with a colony PCR followed by a 
hybridization with (AG) 15 oligonucleotides as probes. This most 
likely helped us to exclude the majority of non-microsateUite 
containing colonies and increased the efTiciency of isolating 
microsatellite-containing clones. Of the 108 clones sequenced, 
97 clones contained microsatellites and only nine (8.3%) did not. 
The sequences of two clones were of poor quality. 

High allelic polymorphism and H^/PIC and simple single-locus 
inheritance patterns of the microsatellites developed in this study 
suggest that the developed markers could be used for various 
genetics, genomics, breeding, DNA fingerprinting, tree forensic 
and genetic resource conservation studies and applications in black 
spruce and red spruce. Indeed six (RPMSA04, RPMSA09, 
RPMSA12, RPMSA13, RPMSA17, and RPMSA33) of the 14 
polymorphic microsatellite markers developed have already been 
mapped to six linkage groups in black spruce [6] . Thus, our study 
provides a rich resource of highly informative single-copy 
microsatellite DNA markers in spruce. 



PLOS ONE I www.plosone.org 



9 



August 2014 | Volunne 9 | Issue 8 | el 03789 



AFLP-Based Microsatellite Markers in Spruce 



S «l E 
-a.S 



Ml 



O D- 

— ni 

u 3t 



o 
u 



<U r-, 



.o ^ 

o3 D 



ii, -I- <u 



.y >, .y 
§ ™ I 

2 

tu — 
Dl Q_ 

^ < 
< 

lii Q 



S, >- Si 

■D -D 

^ — ^ 

LU Q LU 



1:1 ra 



o3 £- 



£ O OJ 

Si Si c 

"D 1^ 3 
QJ O C 

-5^1 



5; S; 



PLOS ONE I www.plosone.org 



10 



August 2014 | Volume 9 | Issue 8 | el 03789 



AFLP-Based Microsatellite Markers in Spruce 



Microsatellite Diversity and Divergence in Black Spruce 
and Red Spruce 

Black spruce individuals had 150 alleles at 16 microsatellite loci 
and 133 alleles at 15 loci after excluding RPMSA05 that showed 
multilocus patterns in red spruce (Table 2). Red spruce individuals 
had 93 alleles at the same 15 SSR loci. Also, the expected and 
observed heterozygosity values were lower in red spruce individ- 
uals. This is consistent with previous reports based on aUozyme 
and cDNA sequence-tagged sites (STS) markers that red spruce 
harbours relatively lower genetic diversity [43-45] as compared to 
sympatric black spruce and white spruce. Genetic diversity at the 
studied microsateUite loci could be regarded as high in black 
spruce and moderate in red spruce accessions studied. However, 
although for a limited number of samples, the microsatellite allelic 
variability results suggest that red spruce cannot be considered as 
genetically depauperate as reported by Perron et al. [44] . 

Red spruce individuals had a number of microsateUite alleles 
(30) that were not detected in black spruce individuals (Table 2; 
Table S3; Table S4). Also locus RPMSA06 was monomorphic for 
one allele in black spruce but had three alleles in red spruce 
individuals. Based on the STS markers. Perron et al. (2000) 
concluded that genetic diversity of red spruce is a subset of black 
spruce genetic diversity [44] . Our microsatellite results, albeit for a 
small but equal number of samples for black and red spruce, do 
not support this conclusion. The microsatellite locus RPMSAOl 
was monomorphic with species-specific alleles in black spruce 
(allele 206) and red spruce (allele 170). If these results hold true for 
a larger number of black spruce and red spruce individuals, the 
RPMSAOl locus could be used to unambiguously dilferentiate 
between closely related black and red spruce. These two spruce 
species hybridize naturally where their ranges overlap, and are 
difficult to differentiate morphologically. Also, no species-specific 
allozyme, ribosomal, mitochondrial or chloroplast DNA markers 
could be found for these species [46-48]. However, RAPD 
markers specific to black and red spruce have been reported [49] . 
But RAPDs are not reliable markers until com erted to codom- 
inant and reliable SCAR markers. Therefore, codominant 
markers, such as the microsatellite DNA markers, could provide 
excellent tools to dilferentiate between black spruce and red 
spruce. 

Inheritance, Segregation and Linkage of Microsatellite 
Markers 

The Mendelian inheritance and segregation of microsatellite 
variants in the progenies of two controlled crosses (Table 3) 
demonstrate that the variants at each of the 1 0 microsatellite loci 
are controlled by a single nuclear locus and that the microsatellite 
variants (alleles) at these loci are codominandy inherited in black 
spruce. Segregation distortions were observed at two microsatellite 
loci: RPMSA04, and RPMSA17 (Table 3). However, for 
RPMSA04, segregation distortion was observed in progeny of 
one controlled cross and not in that of the other. Segregation 
distortion or transmission ratio distortion of molecular markers has 
been reported in many plant species, including spruce [6,50,51]. 
In a black spruce-red spruce BCi backcross, 47% of SAMPL 
(selectively amplified microsatellite polymorphic loci) markers 
showed distorted segregation ratios [50]. In Norway spruce, three 
(2.3%) out of the 133 tests with 51 SSR markers showed a 
significant deviation from MendeUan segregation ratios [51]. In 
other species, for example Isis, approximately 1/3 of all 
microsatellite markers in each linkage map revealed significant 
transmission ratio distortion [52]. In our study, two of 15 (13.3%) 
tests significandy departed from the Mendehan ratios. This level 



was comparable and falls within the range of distortion values 
observed in intraspecific crosses [53]. 

Segregation/ transmission distortion is a common phenomenon 
in plants [54] , with backcrosses and F2 showing high percentage of 
loci showing transmission ratio distortion [54—56]. Several 
biological mechanisms and experimental errors in scoring the 
mapping population genotypes can cause deviations from expected 
segregation ratios. Although we thoroughly verified marker 
scoring, one of the two parents in one of the two cases showing 
segregation distortions had null alleles. This may have caused 
some ambiguity in scoring the progeny genotypes. The presence of 
nuU alleles can cause segregation distortion in the crosses [57]. 
Biological mechanisms that can cause segregation distortion 
include chromosome/chromatin loss [54], divergence of the 
parental genotypes [53,54], inbreeding depression, presence of 
the transmission distorter loci or genomic regions [55], meiotic 
drive locus [56], and genetic load and recessive lethal alleles [58]. 
We examined the inheritance of microsateUite markers in f 2 
progeny of black spruce. The parents of this progeny were related 
by descent. Back spruce shows high inbreeding depression [59]. 
Therefore, the observed segregation distortion may have been 
caused by inbreeding depression. 

Our study in combination with the results reported by Kang 
et al. [6] suggests that there is no linkage among RPMSA04, 
RPMSA07, RPMSA09, RPMSAll, RPMSA12, RPMSA13, 
RPMSA17, RPMSA26, and RPMSA33, and between RPMSA13 
and RPMSA22. Therefore, these microsatellite loci could be 
simultaneously used in population genetic stucUes in black spruce. 

Cross-species Transferability of Microsatellite Markers 

Primers for 14 out of 16 microsatellite loci developed in black 
spruce were able to amplify microsateUites in at least four other 
species, with primers for 1 2 loci yielcUng amplification products in 
aU six spruce species. In red spruce, primers for aU 16 
microsateUite loci worked, resolving single-locus variants for 15 
loci and multUocus patterns for one microsateUite locus (Table 2). 
Thus, it could be inferred that most of the microsateUite markers 
developed from black spruce could potentially be used for various 
genetics, genomics and breeding studies in five other spruce 
species tested. The transferability rate in this study was higher than 
previously reported for spruce [2,15]. It has been suggested that if 
a particular microsateUite locus could be resolved in black, white, 
red, and Sitka spruce by using the same primers, it is likely that the 
primers will be widely transferable to other spruce species as well 
[15]. The high cross-species transferabUity observed in the present 
study indicates that the regions flanking the microsateUite 
sequences are weU conserved in the six spruce species studied. 
The cross species transferability of microsatellite markers in Picea 
species wUl reduce the cost of marker development. 

The success of cross-species transferability of microsatellite loci 
depends on the evolutionary relatedness of the taxa being 
examined. In this study, primers developed for aU 16 microsateUite 
loci from black spruce sequences resolved microsateUites in red 
spruce. This is consistent with close genetic and phylogenetic 
relationships between these two species [44,60]. It is worth noting 
that primers for RPMSA05 resolved a single highly polymorphic 
microsateUite kx us in black spruce but yielded multilocus patterns 
in red spruce, and other four spruce species tested. 

Conclusions 

We have developed a reliable and highly efficient approach for 
isolating and developing highly informative microsateUite DNA 
markers with simple single-locus patterns in black spruce and red 
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spruce. The markers inherited in a Mendelian single-locus 
codominant fashion and were not hnked. Most of the microsat- 
ellite markers could be successfully transferred to Norway, Sitka, 
white, and Engelmann spruce. The new set of microsatellite 
markers developed here could be used for various genetics, 
genomics, breeding, DNA fingerprinting, tree forensics and 
genetic resource conservation studies and applications in six 
spruce and potentially other spruce species. The microsatellite 
markers reported here are based on developing primers for only 34 
out of 94 unique SSR-containing sequences. Thus, additional 
markers could be developed from 60 microsateUite-containing 
sequences that we did not use. 
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Table SI MicrosateUite-containing clones, repeat mo- 
tifs and insert size in black spruce. 

(DOCX) 

Table S2 Microsatellite markers screened; their repeat 
motifs and primer sequences. 

(DOCX) 

Table S3 Allele frequencies at characterized microsat- 
ellite loci in 30 individuals of black spruce {Picea 
mariana) . 

(DOCX) 

Table S4 Allele frequencies at characterized microsat- 
ellite loci in 30 individuals of red spruce {Picea rubens). 

(DOCX) 

Table S5 Joint two-locus segregation patterns and chi- 
square analyses for testing of linkage between micro- 

References 

1. Rajora OP, Rahman MH, Buchert GP, Dancik BP (2000) Micro.satrilitc DNA 
analysis of genetic effects of harvesting in old-growth eastern white pine [Pinus 
strobus) in Ontario. Mol Ecol 9: 339-348. 

2. Rajora OP, Rahman MH, Dayanandan S, Mosseler A (2001) Isolation, 
characterization, inheritance and hnkage of microsatellite DNA markers in white 
spruce [Firea glaucd) and their usefulness in other spruce species. Mol Gen 
Genet 264: 871-882. 

3. Rahman MH, Rajora OP (2001) Microsatellite DNA analysis of somaelonal 
variation in tissue culture micropropagated trembling aspen plants. Plant Cell 
Rep 20: .^Sl-.^Se. 

4. Rahman MH, Rajora OP (2002). Microsatellite DNA fingerprinting, differen- 
tiation and genetic relationships of clones, cultivars and varieties of six poplar 
species from three sections of the genus Populus. Genome 45: 1083-1094. 

5. Rajora OP, Rahman MH (2003) Microsatellite DNA and RAPD fmgerprinting, 
identification and genetic relationships of hybrid poplar {Populus xcanadensis) 
cultivars. Theor Appl Genet 106: 470—477. 

6. Kang BY, Mann IK, Major JE, Rajora OP (2010) Near-saturated and complete 
genetic linkage map of black spruce IPicm mariana). BMC Genomics 11:515. 

7. Fageria MS, Rajora OP (2013) Effects of har\'esting of increasing intensities on 
genetic diversity and population structure of white spruce. Evolutionary 
Applications 6:778-794. 

8. Bashalkhanov S, Eckert AJ, Rajora OP (2013) Genetic signatures of natural 
selection in response to air pollution in red spruce [Picea rubens, Pinaceae). Mol 
Ecol 22: 5877-5889. 

9. Hamblin MT, Warburton ML, Buckler ES (2007) Empirical comparison of 
simple sequence repeats and single nucleotide polymorphisms in assessment of 
maize diversity and relatedness. PLoS One 2: el 367. 

10. Zane L, Bargelloni L, Patamello T (2002) Strategies for microsateUite isolation: a 
review. Mol Ecol 11:1-16. 

11. Cho YG, Ishii T, Temnykh S, Chen X, Lipovich L, et al. (2000) Diversity of 
microsatellites deri\'ed from genomic libraries and GenBank sequences in rice 
(Oryza sativa L.). Theor ,\ppl {,enet 100:713-722. 

12. Thiel T, Miehalek W, Varshney RK, Graner A (2003) Exploiting EST databases 
for the de\ elopment and characterization of gene-derived SSR-markers in barley 
{Hordeum vulgare L.). Theor Appl (ienet 106: 411-422. 

13. Scott KD, Eggler P, Seaton G, Rossetto .M, Ablett EM, et al. (2000) Analysis of 
SSRs derived from grape ESTs. Theor Appl Genet 100:723-726. 



satellite DNA loci in the 46x14 cross of black spruce 
{Picea mariana). 

(DOCX) 

Table S6 MicrosateUite-containing sequences of black 
spruce used for primer design. 

(DOCX) 

Acl(nowledgments 

Wc thank John Major for providing the black spruce pedigreed material 
that was originally produced by late Dr. Kris Morgenstern and Dr. Tim 
Boyle, Alex Mosseler lor prox'iding red spruce samples, Michael Stoehr for 
providing .Sitka spruce samples, Howard Frame for providing Norway 
spruce samples, and Joy Roy, AUan Bland and Nicholas Buckley for 
laboratory assistance with microsatellite genotyping. 

Author Contributions 

Conceived and designed the experiments: YS OPR. Performed the 
experiments: YS NF. Analyzed the data: YS NF OPR. Contributed 
reagents/ materials/ analysis tools: OPR. Wrote the paper: OPR YS NF. 
Contributed to the conception of idea for developing microsatellite markers 
from microsatellite-enriched AFLP fragments, developed and optimized 
the protocol, developed and optimized microsatellite markers and 
contributed to the manuscript writing and revision: YS. Contributed to 
the lab work and data scoring and analysis for determining the 
informtaiveness, inheritance and cross-species transferability of the 
microsatellite markers, and to the manuscript preparation: NF. Principal 
Investigator of the project and contributed to the study conception, 
experimental design, data analysis, overall research directions and 
guidance, and manuscript writing and revision, and obtained funding 
and provided laboratory facilities: OPR. 



14. Liewlaksaneeyanawin C, Ritland CE, El-Kassaby YA, Ritland K (2004) Single- 
copy, species-transferable microsatellite markers developed from loblolly pine 
ESTs. Theor Appl Genet 109: 361-369. 

15. Rungis D, Berube Y, ZhangJ, Ralph S, Ritiand CE, et al. (2004) Robust simple 
sequence repeat markers for spruce {Picea spp.) from expressed sequence tags. 
Theor Appl Genet 109: 1283-1294. 

16. Murray BG (1998) Nuclear DNA amounts in g^Tnnosperms. Ann Bot 82 (suppl 
1): ,3-15. 

17. Kinlaw GS, Nealc DB (1997) Complex gene families in pine genomes. Trends 
Plant Sci 2: 356-359. 

18. Nystedt B, Street NR, Wetterbom A, Zuccolo A, Lin YC, et al. (2013) The 
Norway spruce genome sequence and conifer genome evolution. Nature 
497(745 1):579-584. 

19. Birol I, Raymond A, Jackman SD, Pleasancc S, Coope R, et al. (2013) 
Assembling the 20 (}b white spruce [F/rea glauca) genome from whole-genome 
shotgun sequencing data. Bioiitlormaties 29:1492-1497. 

20. Pfciller A, 01i\ ieri .Wl. Morgante M (1997) Identification aitd characterization 
of microsatellites in Xorwav spruce [Ficea ahies K.). (ienome 40: 411-419. 

21. Hodgetts RB, Aleksiuk MA, Brown A, Clarke C, Maedonald E, et al. (2001) 
Development of microsatellite markers for white spruce {Picea glauca) and 
related species. Theor Appl Genet 102: 1252-1258. 

22. Scotti I, Magni F, Paglia GP, Morgante M (2002) Trinucleotide microsatellites 
in Norway spruce {Picea abies): their features and the development of molecular 
markers. Theor Appl Genet 106: 40-50. 

23. Li Y, Liang L, Ge XJ (2010) Development of microsatellite loci for Pinus 
komifnsis (Pinaceae). Am J Bot 97: e39-c41. 

24. Georgi LL, Wang Y, Yvergniaux D, Ormsbce T, Inigo M, et al. (2002) 
Construction of a BAG library and its application to the identification of simple 
sequence repeats in peach [Frunus persica (L.) BatschJ. Theor Appl (jcnet 105: 
1151-1158. 

25. Chen H, Pulido JC, Duyk CM (1995) MATS: a rapid and efficient method for 
the development of microsatellite markers from YACs. Genomics 25: 1—8. 

26. Elsik CG, Wilhams CG (2001) Low-copy microsatellite recovery from a conifer 
genome. Theor Appl Genet 103: 1189-1195. 

27. Zhou Y, Bui T, Auckland LD, Williams CG (2002) Undermethylated DNA as a 
source of microsatellites from a conifer genome. Genome 45: 91—99. 

28. Rajora OP, Pluhar SA (2003) (ienetie diversity impacts of forest fires, forest 
harvesting and alternative reforestation practices in black spruce {Picea 
mariana). Theor Appl Genet 106: 1203-1212. 



PLOS ONE I www.plosone.org 



12 



August 2014 I Volume 9 | Issue 8 | e103789 



AFLP-Based Microsatellite Markers in Spruce 



29. Yamamoto T, Kimura T, Shoda M, Ban Y, Hayashi T, ct al. (2002) 
Development of microsatellite markers in the Japanese pear [Pyrus pyrifolia 
Nakai). Molecular Ecology Notes 2: 14-16. 

30. Yamamoto T, Mochida K, Imai T, Shi YZ, Ogiwara I, et al. (2002) 
Microsatellite markers in peach [Prunus persica (L.) Batsch] derived from an 
enriched genomic and cDNA libraries. Molecular Ecology Notes 2: 298-301. 

31. Rozcn S, Skalctsky H (2000) Primer 3 on the WWW for general users and for 
hioldiiisl programmers. In: Clifton XJ. editor. Methods in Molecular Biology. 
Humana Press, Totowa, NJ: pp 365 386. 

32. Cue LM, iMace ES, Crouch JH, Ouang YD. Long TD, rt al. (2008) Isolation 
and characterization ol novel microsatellite markers and their application for 
diversity assessment in cultivated groundnut [Arachis hypogaea). BMC Plemt Biol 
8: 55. 

33. Sokal RR, Rohlf FJ (1995) Biometr\'. Third edition. W.H. Freeman and 
Company, New York. 

34. Dobrezeniecka S, Xkongolo KK, Michael P, Wyss S, Mehcs M (2009) 
Identification and charartcrizaiiou ol' microsatellite markers for genetic analysis 
of black spruce [Picea mariana (Mill.) populations. Silvae (ienetica 58:168-172. 

35. Taramino G, Tarchini R, Ferrario S, Lee M, Pe' ME (1997) Characterization 
and mapping of simple sequence repeats (SSRs) in Sorghum bicolor. Theor Appl 
Genet 95: 66-72. 

36. Areshchenkova T, Ganal MW (2002) Comparative analysis of polymorphism 
and chromosomal location of tomato microsatellite markers isolated from 
different sources. Theor Appl Genet 104: 229-235. 

37. Donini P, Stephenson P, Bryan GJ, Koebner RMD (1998) The potential of 
microsatellites for high throughput genetic diversity assessment in wheat and 
barley, (ienetie Resources and Crop Evolution 45: 415—421. 

38. He GH, Meng RH, Newman M, Gao GQ, Pittman RN, et al. (2003) 
Microsatellites as DNA markers in cultivated peanut {Arachis hypogaea L.) BMC 
Hant Biol 3:3. 

39. Scotti I, Paglia GP, Magni F, Morgante M (2002) Efficient development of 

dinucleotide microsatellite markers in Nor\vay spruce {Picea abies Karst.) 
through dot-blot selection. Theor Appl Genet 104: 1035-1041. 

40. Rabinowicz PD, Citek R, Budiman MA, Andrew Nunberg A, Bedell JA, et al. 
(2005) Differential mcthylation of genes and repeats in land plants. Genome Res 
2005 15: 1431-1440. 

41 . YangJB, Li HT, Li DZ, Liu J, Gao LM (2009) Isolation and characterization of 
microsatellite markers in the endangered species Taxus wallichiana using the 
FIASCO mediod. HortScience 44:2043-2045. 

42. Miao YC, Lang XD, Li SF, Su JR, Wang YH (2012) Characterization of 15 
polvmorphic microsatellite loci for Cephalotaxus oliveri (Cephalotaxaceae), a 
conifer of medicinal importance. IntJ Mol Sci 13:1 1165-11 172. 

43. Hawlev (ij, DeHayes DH (1994) Cjcnetic di\'ersity and population structure of 
red spruce {Picea ruhens). CanJ Bot 72: 1778-1786. 

44. Perron M, Perry DJ, Andalo C, Bousquet J (2000) Evidence from sequence- 
tagged-site markers of a recent progenitor-derivative species pair in conifers. 
Proc Natl Acad Sci USA 97: 11331-11336. 

45. Rajora OP, Mosseler A, Major JE (2000) Indicators of population viability in red 
spruce, Picea ruhens. II. Genetic diversity, population structure, and mating 
behaviour. CanJ Bot 78: 941-956. 



46. Eekert RT (1989) Genetic variation in red spruce and its relation to forest 
decline in northeastern United States. In: Bueher JB, Bucher-Wallin I (eds) Air 
pollution and forest decline. lUFRO, Birmensodorf, pp 319-324. 

47. Bobola MS, Eckert RT, Klein AS (1992) Restriction fragment variation in the 
nuclear ribosomal DNA repeat unit within and between Picea rubens and Picea 
mariana. CanJ For Res 22: 255-263. 

48. Bobola MS, Eckert RT, Klein AS, Strapelfeldt K, Smith DE, et al. (1996) Using 
nuclear, and organelle DNA markers to discriminate among Picea rubens, Picea 
mariana, and their hybrids. CanJ For Res 26: 433^43. 

49. Perron M, Gordon AG, Bousquet J (1995) Species-specific RAPD fingerprints 
for the closely related Picea mariana and Picea ruhens. Theor Appl (ienet 91: 
142-149. 

50. Yazdani R, Scotti l.Jansson Cj, Plomion C, Mathur G (2003) Inheritance and 
diversity of simple sequence repeat (SSR) microsatellite markers in various 
families of Picea abies. Hereditas 138: 219-227. 

51. Kang BY, MajorJE, Rajora OP (201 1) A high-density genetic linkage map of a 
black spruce {Picea mariana) xrcA spruce {Picea rubens) interspecific hybrid. 
Genome 54:128-143. 

52. Tang SX, Okashah RA, Knapp SJ, Arnold ML, Martin NH (2010) 
Transmission ratio distortion results in asymmetric introgression in Louisiana 
Ins. BMC Plant Biol 10:48. 

53. Jenezewski E, (iherardi M, Bonnin I, ProsperiJM, Olivieri I, et al. (1997) Insight 
on segregation distorlions in two intraspeeilie crosses between annual species of 
Medicago (Leguminosae). Theor Appl Genet 94: 682-691. 

54. Zamir D, Tadmor Y (1986) Unequal segregation of nuclcEir genes in plants. Bot 
Gaz 147: 355-358. 

55. Fishman L, Aagaard J, Tuthill JC (2008) Towards the evolutionary genomics of 
gametophytic divergence: patterns of transmission ratio distortion in monkey 
flower {Mimulus) hybrids reveal a complex genetic basis for conspeciflc pollen 
precedence. Evolution 62: 2958-2970. 

56. Fishman L, Willis JH (2005) A novel meiotic drive locus almost completely 
distorts segregation in Mimulus (monkeyflower) hybrids. Genetics 169: 347-353. 

57. Reece KS, Ribeiro WL, Gaffney PM, Carnegie RB, Allen SK Jr (2004) 
Microsatellite marker development and analysis in the eastern oyster [Crassostrea 
virginica): confirmation of null alleles and non-Mendelian segregation ratios. 
J Hered 95; 346-352. 

58. Bradshaw HD Jr, Stettler RF (1994) Molecular genetics of growth and 
development in Populus. 11. Segregation distortion due to genetic load. Theor 
Appl Genet 89: 551-558. 

59. Morgenstern EK (1972) Preliminary estimates of inbreeding in natural 
populations of black spruce, Picea mariana. CanJ Genetics Cytol 14:443-446. 

60. Ran JH, Wei XX, Wang XQ_(2006) Molecular phylogeny and biogeography of 
Picea (Pinaceae): Implications for phylogeographcal studies using cytoplasmic 
haplotypes. Mol Phylogenet Evol 41:405-419. 

61. Scotti I, Magni F, Fink R, Powell W, BmeUi G, et al. (2000) MierosateUite 
repeats are not randomly distributed within Norway spruce [Picea abies K.) 
expressed sequences. Genome 43: 41—46. 



PLCS ONE I www.plosone.org 



13 



August 2014 I Volume 9 | Issue 8 | e103789 



